Nonlinearity induced destruction of resonant tunneling in the Wannier-Stark problem 
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We present detailed numerical results on the dynamics of a Bose-Einstein condensate in a tilted 
periodic optical lattice over many Bloch periods. We show that an increasing atom-atom interaction 
systematically affects coherent tunneling, and eventually destroys the resonant tunneling peaks. 

PACS numbers: 03.75.Lm,03.65.Xp,05.60.Gg 



Experiments with cold and ultracold atoms made it 
possible in the last decade to prepare and control the 
centre-of-mass motion of atoms with unprecedented pre- 
cision. Many toy models of either many-body solid state 
physics Q, 0, 0, 0, H or of simple Hamiltonian systems, 
whose complexity arises from an external driving force 
0, were realized with the exceptional control offered by 
static or time-dependent optical potentials. 

Particularly Bose-Einstein condensates (BEC) whose 
initial momentum spread can be adjusted in width and 
absolute position have proved to be an extremely helpful 
experimental tool USE 110 In addition, a BEC offers 
interesting new features originating from the intrinsic in- 
teractions between the atoms. Examples of such effects 
are new quantum phases |Tcj . soliton-like motion 
the occurrence of energetic or dynamical instabilities in 
condensates Q, 03> EH > ° r the decay and subsequent 
revival of Bloch oscillations (BO) fidj] . 

We focus on the evolution of a BEC loaded into a one- 
dimensional lattice and subjected to an additional static 
force F, which is most easily realized and controlled by 
accelerating the optical lattice U Q- In previous ex- 
periments, a BEC was accelerated to allow for a single 
crossing of the Brillouin zone (BZ), and two effects were 
observed: for large accelerations, an enhanced tunnel- 
ing probability from the ground state band to the first 
excited band due to the atom-atom interaction was mea- 
sured 0, ll| • Secondly, for smaller accelerations (where 
tunneling is negligible) signatures of a dynamical insta- 
bility in the BEC were observed By contrast, here 
we investigate the dynamics of a BEC performing many 
Bloch oscillations (BO), and we ask ourselves how the 
atom-atom interaction affects tunneling for a sequence 
of BZ crossings. In particular, we scan F to study the 
impact of the atom-atom interaction on resonantly en- 
hanced tunneling (RET), for which the standard Landau- 
Zener prediction is modified even in the absence of inter- 
actions 0| . The RET leads to a faster decay of the Wan- 
nier states trapped in the potential wells. With the sur- 
vival probability and the recurrence probability (see Eqs. 
Q and © below) we present two consistent measures for 
the nonlinear RET which define clear experimental signa- 
tures of the destruction of the coherent tunneling process 
inside the periodic potential. 

If we neglect interactions for a moment, our system 



will be described by the Hamiltonian 
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Here is the spatial period of the optical lattice with 
maximal amplitude V, and M the atomic mass. Eq. l|T| 
defines the well-known Wannier-Stark problem, which 
gives rise to BO with period Teioch = h/di,F (h is 
Planck's constant). If tunneling is small, we can view 
the system as moving at a constant speed in momentum 
space within the fundamental BZ. At the zone edge, most 
of the wave packet is reflected (giving rise to BO) while a 
small part can tunnel across the first band gap to the next 
higher-lying energy band and then escape quickly by suc- 
cessive tunneling events across the smaller (higher) band 
gaps. Landau-Zener theory predicts a decay rate [la] 



T(F) oc Fe- 
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where b is proportional to the square of the energy gaps. 
Eq. (J2J is modified by RET which occurs when two 
Wannier-Stark levels in neighboring potential wells are 
coupled strongly due to their accidental degeneracy. The 
RET results in pronounced peaks in the tunneling rates, 
e.g., as a function of 1/F, on top of the global exponential 
decay described by 16] . In this paper we investigate 
the impact of the effective shift of the Wannier-Stark lev- 
els by a nonlinear interaction term. 

For the linear problem l(T|l. the decay rates have been 
measured previously in the regime of short life-times in 
the ground state band (of the order of 100 /is), where 
T(F) is essentially smooth ^]J. Since RET is a coherent 
quantum effect, the peaks should be sensitively affected 
by the atom-atom interaction, which can be varied exper- 
imentally by changing either the density of the BEC or 
through the atom-atom scattering potential via a Fesh- 
bach resonance 0|. Our results are a consequence of 
many sequential Landau-Zener events, and they show the 
destruction of a RET peak with increasing interaction 
strength, in a regime which is experimentally accessible. 

We use a fully 3D Gross-Pitaevskii equation (CPE) 
[l9| to describe the temporal evolution of a BEC which 
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is subject to realistic potentials: 
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ip(r,t) represents the condensate wave function, and the 
frequencies oj x and uj r characterize the longitudinal and 
transverse harmonic confinement (here with cylindrical 
symmetry: p — \/ y 2 + z 2 ). We fixed — 1-56 /im and 
V/E-r = 5 for our computations, with the recoil energy 
En = p 2 l /2M for pr = hn/d,L, and the recoil period 
Tr = h/E-R. The above values for dr and V were re- 
alized in the experiments reported in 0, 0, [t| based on 
two laser beams propagating at an angle different from 
7T. In Eq. @ , the nonlinear coupling constant is given by 
g = Airh 2 a s /M , where a s is the s-wave scattering length 
and N the number of atoms in the BEC 0, |2(| . The 
dimensionless nonlinearity C — <?71o/(8£r) is computed 
from the peak density of the initial state of the conden- 
sate, with C — 0.027... 0.31 for the experimentally in- 
vestigated range of 0, and with C = 0.5 reached in [2l| . 
Here we focus on C > 0, but report briefly also on at- 
tractive interactions with C < 0. The latter case leads to 
a fundamentally different behavior of the system because 
the collapse of the condensate introduces an additional 
time scale, which for experimentally relevant parameters 
is of the order of 10 msec jl8l |22| (slightly longer than 
^Bioch = 1.8 ... 3.0 msec here). 

The GPE © is numerically integrated using finite dif- 
ference propagation, adapted by a predictor-corrector es- 
timate to reliably evaluate the nonlinear interaction |19|| . 
Since our system is essentially the problem of a constantly 
accelerated particle for the part of the wave function 
which has tunnelled out of the first BZ already, one must 
be careful with the application of absorbing boundary 
conditions or complex coordinate methods |23l l24j . To 
avoid any spurious effects due to the fast spreading, we 
use a large numerical basis. In this way, we fully cover 
the 3D expansion of the entire wave packet, including its 
tunnelled tail, without the use of non-Hcrmitian poten- 
tials. The initial state propagated by J3J is the relaxed 
condensate wave function, adiabatically loaded into the 
confining potential given by the harmonic trap and the 
optical lattice (with F = 0). Approximate analytic forms 
of the relaxed state are found, e.g., in [25|, but we used 
an imaginary time propagation to reliably compute the 
initial state for C > 0. 

The linear decay rates for non-interacting atoms in the 
optical lattice are computed from the spectrum of the ID 
Wannier-Stark problem of Eq. Q using, e.g., the method 
of nu. Those linear rates are plotted in Fig. ^ The 
maxima in the rates occur when Fd^m (with m integer) 
is close to the difference between the first two energy 
bands (averaged over the BZ) of the F — problem [lq . 
The actual peaks are slightly shifted with respect to the 
above estimate (marked by arrows in the inset of Fig.^l, 
owing to a field-induced level shift pl|. 
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FIG. 1: (color online). Tunneling rates obtained by expo- 
nential fits to data of P SU r(i) as the solid-line fits in Fig. [3] 
Here the peak in the box in the inset is scanned locally, while 
globally the rates follow an exponential law (dashed line in 
the inset). C = —0.31 (diamonds), C — —0.065 (crosses), 
C = (solid line), C = 0.027 (pyramids), C = 0.065 (stars), 
C — 0.12 (circles), C = 0.31 (squares). The "error" bars 
interpolate different exponential fits such as the dot-dashed 
ones in Fig. |S| The arrows in the inset mark the peak posi- 
tions as predicted by the simple argument stated in the text. 



Experimentally, the most easily measurable quantity is 
the momentum distribution of the BEC obtained from a 
free expansion after the evolution inside the lattice. From 
the momentum distribution we determine the survival 
probability by projection of the evolved state ip{Pj t) onto 
the support of the initial state 



dp x 



dp y dp z \ip{p,t)\- 
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where p c > 3pjt is a good choice since three momentum 
peaks are initially significantly populated, corresponding 

to-2p R ,o,2pR Sea- 

Figure [21 shows the initial population in momentum 
space (inset in (a)) as compared with the population after 
10 BO periods, for both the linear and the nonlinear 
case. The increase of C > has two effects: firstly, it 
enhances the tunneling for the first few crossings of the 
BZ. Secondly, it scrambles the out-coupled part of the 
wave function (see Fig. [2 and its complement in Fig. 0] 
below), as previously observed in (3- 00]- The change in 
the momentum distributions after various Landau- Zener 
events is a manifestation of the intrinsic instability of the 
nonlinear GPE dynamics |9Lll2T|. 

Instead of studying the details of the distributions 
shown in Fig. |3 we will focus on the temporal decay 
of the survival probability in the following. Figure |3| 
presents P S m(t), which for the linear case has an expo- 
nential form (apart from the t — > limit 26] ) 



(5) 



3 




FIG. 2: Momentum distributions after 10 BOs, for 
TR/TBioch = 1.428. C = (dotted) compared with (a) 
C — 0.027 (full line; the inset shows the corresponding t = 
distribution), (b) C = 0.065, (c) C = 0.12, and (d) C = 0.31. 



with the characteristic exponent Y . The temporal behav- 
ior of P aux depends significantly on C. For C — ±0.31, 
we observe clear deviations from a purely exponential de- 
cay, as present for small C . A repulsive nonlinearity ini- 
tially enhances the tunneling more than after about five 
crossings of the BZ (see fits to data in Fig. [3J. This de- 
viation from the mono-exponential behavior means that 
the tunneling events occurring at different integer mul- 
tiples of the Bloch period are correlated by the presence 
of the nonlinearity. Since the remaining density becomes 
smaller, the impact of the nonlinearity becomes less. The 
result is that the rate Y is defined only locally in time, 
and its value systematically decreases as time increases. 

An attractive interaction can stabilize the system at the 
RET peak, which is shown for C = -0.31 in Fig. 3(b). 
For optimal comparison, we chose the same initial state 
(for C = +0.31) which then was evolved for F ^ with 
C = —0.31. Such a scenario could be realized by a sud- 
den change of the sign of the scattering length through 
a Feshbach resonance 0] . This result is consistent with 
studies of simpler models, where a resonance state can be 
stabilized at system-specific strengths of the nonlinearity 
IH|23- 

The impact of the nonlinearity on the dynamical evolu- 
tion of the "closed" system confined to the fundamental 
BZ can be studied with the help of the recurrence prob- 
ability i defined by the autocorrelation 

P rec (i) = |(^(t)|V>(t = 0)>| 2 . (6) 

The BO manifest themselves as the periodic oscillations 
in .P r ec(i) plotted in Fig. These oscillations are less and 
less pronounced with increasing C, in much the same way 
as the momentum peaks are washed out when the band 
edge is crossed in the regime of instability 9] . In contrast 
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FIG. 3: (color online). P Bur (t) for (a) r R /T B ioch = 1-428 
and (b) 1.613 (the peak maximum in Fig. 0. C is scanned 
from -0.31 (diamonds in (b) only), (solid line), 0.027 (dot- 
ted line), 0.065 (dashed line), 0.12 (circles) to 0.31 (squares). 
The gray/red solid lines show global exponential fits to the 
C = 0.31 data, while the dot-dashed lines show exponential 
fits for small and large t respectively. From those fits, the 
rates in Fig. Q and their systematical variation in time are 
obtained. The step-like structures reflect the periodic BO 
and are correlated with the dephased oscillations in Fig. |1] 

to the survival probability, P mc is a phase sensitive mea- 
sure, and therefore it shows - in addition to the temporal 
decay - the dephasing of the BO due to the nonlinearity. 
For C — 0, the recurrence maxima decay in time with the 
same rate as -P S ur(i), which offers an alternative method 
for extracting Y. For C ^ 0, P Tec can be integrated over 
time, and the rates are extractable by the approximate 
proportionality between the integrated area and the in- 
verse decay rate (recalling that J dtf (t) exp(-tY) ~ 1/Y 
to leading order, for a periodic function f(t)). The latter 
approach works because we can determine the linear rate 
from a direct fit to P Tec and then compare the ratio of the 
linear and the nonlinear area (denoted by Aq and Ac). 
This rough estimate Yc ~ YqAq/Ac agrees within 25% 
with the rate extracted from the fits to the data of Fig. [3] 
The estimate could be improved if we knew the analytic 
form of the function f(t), and it breaks down for large 
C, when the periodic oscillations in P rec are destroyed. 

Having introduced two methods to extract the tunnel- 
ing rates, we scan the parameter F across a RET peak of 
the globally exponential curve Y(l/F) (see Fig. The 
scanned range in F corresponds to values of lattice ac- 
celerations between 0.99 ms~ 2 and 1.65 ms~ 2 , which are 
standard in experiments 0, || . 

A repulsive nonlinearity particularly affects the wings 
of the peak and, for small C, much less the peak maxi- 
mum. The global increase of Y with increasing C is qual- 
itatively predicted in |28| , with enhanced single Landau- 
Zener crossing probabilities induced by the effective re- 
duction of the energy gap due to the nonlinearity. The 
left and right-most points in Fig. ^ are m the regime 
where an amended version of indeed applies |2g , and 
here Y/F is approximately proportional to C. However, 
near the peak, the rates do not follow a simple scaling 
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FIG. 4: (color online). P rC c(i) for the data shown in Fig. OJb) 
with C = (circles), 0.027 (dashed), 0.12 (thin line), and 0.31 
(thick line). The dot-dashed line presents an exponential fit 
to the maxima of the C = data. 

law as a function of C, and the argumentation of 
does not apply. 

For C < we also observe the destruction of the RET 
peak. For C — —0.065, the BEC clearly stabilizes in 
the potential wells, whilst for C — —0.31 the situation 
is more complicated (see Fig. The precise dynam- 
ics of the system is governed by the two separate time 
scales for tunneling and collapse, which strongly depend 
on parameters in the sensitive RET regime. 

In an experiment, lu x can cither be set to zero or de- 
creased to u> x /2ir < 1 Hz to realize a quasi ID nonlinear 



Wannier-Stark problem. We verified that letting oj x tend 
to zero for the evolution with F 0, or applying a small 
finite u> x gives the same results for the BO cycles studied 
here. Furthermore, for < C < 0.05, using the renormal- 
ized nonlinearity of [2(| we observed that a ID version 
of Eq. © reproduces well the 3D data. If |C| is larger, 
the nonlinearity couples the longitudinal and transverse 
degrees of freedom, which affects the dynamics of a real 
BEC in a non-trivial way 0|. The ID computations are 
feasible up to 100 Bloch periods, and this would allow 
one to extract the tunneling rates more reliably. The 
effect of the nonlinearity is, however, hardly visible for 
< C < 0.05, and quantitative predictions for a broad 
range of C relied on 3D computations. 

To summarize, we observed and quantified the defor- 
mation and destruction of the RET peaks due to interac- 
tions in a BEC in an accelerated optical lattice. Our re- 
sults complement ongoing studies of interaction-induced 
processes such as dynamical instabilities or the decay and 
subsequent revival of BO. In the regime of small nonlin- 
earity, where dynamical instabilities are not fully devel- 
oped, the survival and recurrence probabilities experience 
an exponential decay modified by the condensate nonlin- 
earity. The temporal decay of these observables remains 
a useful indicator also for large nonlinearity, even if the 
resonant structure in the tunneling rate is washed out. 

We thank M. Cristiani and D. Ciampini for helpful dis- 
cussions and the Humboldt Foundation (Feodor-Lyncn 
Program), MIUR COFIN-2004, and ESF (QUDEDIS) 
for support. 



[1] M. BenDahan et al, Phys. Rev. Lett. 76, 4508 (1996); 

S.R. Wilkinson et al, ibid. 76, 4512 (1996). 
[2] B.P. Anderson and M.A. Kasevich, Science 282, 1686 

(1998). 

[3] O. Morsch et al, Phys. Rev. Lett. 87, 140402 (2001); M. 

Cristiani et al, Phys. Rev. A 65, 063612 (2002). 
[4] J. Hecker-Denschlag et al, J. Phys. B 35, 3095 (2002). 
[5] G. Roati et al, Phys. Rev. Lett. 92, 230402 (2004). 
[6] M.G. Raizen, Adv. At. Mol. Opt. Phys. 41, 43 (1999); 

D.A. Steck, W.H. Oskay, and M.G. Raizen, Science 293, 

274 (2001); W.K. Hensinger et al, Nature (London) 412, 

52 (2001). 

[7] L. Fallani et al, Phys. Rev. Lett. 93, 140406 (2004). 
[8] M. Jona-Lasinio et al, Phys. Rev. Lett. 91, 230406 
(2003). 

[9] M. Cristiani et al, Opt. Express 12, 4 (2004). 
[10] M. Greiner et al, Nature (London) 415, 39 (2002). 
[11] B. Eiermann et al, Phys. Rev. Lett. 92, 230401 (2004). 
[12] B. Wu and Q. Niu, New J. Phys. 5, 104 (2003). 
[13] Y. Zheng, M. Kostrun, and J. Javanainen, Phys. Rev. 

Lett. 93, 230401 (2004). 
[14] A.R. Kolovsky, Phys. Rev. Lett. 90, 213002 (2003); A. 

Buchleitner and A.R. Kolovsky, ibid., 91, 253002 (2003); 

Q. Thommen, J.C. Garreau, and V. Zehnle, ibid., 91 

210405 (2003). 



[15] Instability at the band edge has also been studied using 

a different experimental protocol in |2|. 
[16] M. Gliick, A.R. Kolovsky, and H.J. Korsch, Phys. Rev. 

Lett. 83, 891 (1999); Phys. Rep. 366, 103 (2002). 
[17] C.F. Bharucha et al, Phys. Rev. A 55, R857 (1997). 
[18] J.L. Roberts et al, Phys. Rev. Lett. 86, 4211 (2001). 
[19] E. Cerboneschi et al, Phys. Lett. A 249, 495 (1998); S. 

Wimberger et al, Phys. Rev. Lett. 94, 130404 (2005). 
[20] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 
[21] N. Gemelke et al, cond-mat/0504311 
[22] S. Wiister, I.J. Hope, and CM. Savage, Phys. Rev. A 71, 

033604 (2005). 

[23] T. Paul, K. Richter, and P. Schlagheck, Phys. Rev. Lett. 

94, 020404 (2005). 
[24] P. Schlagheck and T. Paul, |cond-m at/0402089 N. Moi- 

seyev and L.S. Cederbaum, cond-mat/0406189 
[25] P. Pedri et al, Phys. Rev. Lett. 87, 220401 (2001). 
[26] S.R. Wilkinson et al, Nature (London) 387, 575 (1997). 
[27] D. Witthaut, S. Mossmann, and H.J. Korsch, J. Phys. A 

38, 1777 (2005). 
[28] D.I. Choi and Q. Niu, Phys. Rev. Lett. 82, 2022 (1999); 

O. Zobay, B.M. Garraway, Phys. Rev. A 61, 033603 

(2000). 



